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ABSTRACT 

Relativistic electron degeneracy effects are dominant in ultra-dense plasmas 
(UDP), such as those found in white dwarfs. These effects can be treated sys- 
tematically by obtaining an expansion of the screening length in inverse powers 
of h 2 . In general, our theory leads to an (9(10) effect on the enhancement of 
fusion rates in white dwarfs. Further, it is shown analytically for these stellar 
conditions that Bose statistics of nuclei have a negligible effect on the screening 
length, in consonance with Monte Carlo simulations found in literature. 

Subject headings: screening, fusion, plasma, quantum corrections 
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Introduction 



The evolution of a white dwarf is marked by fusion reactions between carbon and oxy- 
gen nuclei. These reactions are enhanced by the screening of the Coulomb interaction 
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) . It is difficult to gauge the accuracy of these 



calculations. Our goal here is to provide a systematic basis for the estimation of the fusion 
enhancement factor in white dwarfs. 

It is possible that the techniques developed in this pa per may be applicable to studying 

other elementary particle processes in dense stellar interiors (lltoh et al.ll2002l ; lLanganke and Martinez-Pined' 
2003h . 



Nuclei in white dwarfs are submerged primarily in a sea of electrons. This is because the 
probability of finding a nucleus, say, in a two-component UDP (one specie of nuclei of charge 
Ze, and electrons) is down by a factor of 1/Z, compared to finding an electron. Consequently, 
we might expect that Coulomb screening will be dominated by the sea of electrons. These 
electrons are not only degenerate, but also relativistic since the associated fermionic chemical 
potential is much larger than the thermal energy. Physically, the question of the dominant 
mode of screening boils down to how closely these energetic, relativistic electrons surround 
a given nuclei, compared to the proximity of another nucleus. Coulomb repulsion will keep 
charged nuclei approximately a Debye-Huckel length away from each other. We find in this 
paper that the energetic nature of these relativistic electrons keeps them away farther than 
one might expect on naive grounds, approximately a Debye length as well. 



Our earlier paper (IChitanvisll2007l ) was directed towards understanding the scre ening ef- 
fect of electronic quantum fluctuat ions on fusion reactions near the center of the sun ((Salpeter 
1954 ; iGruzinov and Bahcalllll998l ). That paper showed quantum effects are negligible, via a 
systematic expansion of the screening leng th in powers of fi 2 , putting to rest a controversy 
that has ebbed and flowed over the years (lBahcall et al.ll2002l ). In this paper we consider 
a plasma where relativistic quantum effects dominate, which is quite opposite to the solar 
plasma. And so we contemplate the theory in our previous paper in inverse powers of h 2 , 
which will apply when the effects of electron degeneracy dominate. 

This paper also provides an an alytical underpinning to the numerical techniques used 
in the past to study fusion in UDP (jOgatalll997l ; iPollock and Militzer! l2004j ) . These papers 
show that the effect of indistinguishability of nuclei on the screening of fusion rates is small. 
Further, our technique provides an alternative, systematic method for the calculation of the 
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enhancement factors of fusion rates in UDP. Our estimates of enhancement are consistent 



with those given by Ichi maru and Kitamurafllchimaru and Kitamuralll999l ). and with recent 
results of Gasques et a KIGasques et al.ll2005l ). 



2. Screening Formalism for relativistic degenerate electrons 



Let us begin with the classical Poisson-Boltzmann equation for a single species of ions 
and electrons: 



- V 2 = Aixp 

P = P+ + P- 
p + = e n Z exp(—Zecf)/kBT) 

p_ = — e n Z exp(e(p/k B T) (1) 

where e is the magnitude of the electronic charge, k B is Boltzmann's constant, n is the 
average number density, Ze is the ionic charge, and T is the temperature of the system. We 
shall work in the linear regime, by retaining only terms first order in 0, leaving to the next 
section a discussion of nonlinear terms resulting from the Boltzmann distribution: 



^Ann(Z 2 n + Zn)e 2 



k B T 



Ao 




k B T 



47me 2 (Z 2 + Z) 



(2) 



where A is the classical screening length. 



We shall now generalize this method to one where nu clei are treated classically, but elec- 
trons are treated quantum mechanically (jChitanvisl 120071 ) . The quantum-mechanical version 
of the linearized Poisson-Boltzmann equation for a single species of ions and electrons may 
be written in analogy with Eqn. [TJ 
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- V 2 = 4?rp 

P = P+ + P- 

p_ =-e|^({r})| 2 (3) 

where ■?/> is the many-body quantum wave-function for electrons, and {r} refers collectively 
to the electrons in the system, and is the electrostatic potential. 

We now invoke the following scaled variables, in order to ease subsequent calculations: 



A 



knT 



4irZ 2 ne 2 
r 

A 

Ze 2 



r' = (4) 
Ak B T v ; 

where T' is defined differently from the usual plasma parameter. Notice also that the scalar 



potential has been scaled differently than in our previous paper ( IChitanvisI 120071 ). This is 
done to allow for a convenient analysis of higher order contributions in the next section. A 
is the average distance maintained between nuclei. 

Note that the first of EqnsQ] shows that we are using k B T as the energy scale. The 
electrostatic potential is then given by: 

V' 2 0= (0 + 47rr>| 2 -l) (5) 

where if) is a four-vector corresponding to the Dirac spinor. This equation may be obtained 
from a Lagrangian density: 



Co = -^| 2 -v(W) 
v(M) =^ 2 + 4vr0r'H 2 -0 (6) 
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The corresponding Hamiltonian density can be easily derived: 

H = l\V4>\ 2 + v(4>,4>) (7) 

We will now introduce second-quantized notation to deal with the statistics of electrons: 

v(4>, $) - v(4>, j>±) = ^4> 2 -4> + 47r0r'(^) (8) 

where ip are Grassmann variables corresponding to Dirac electrons. The co-existence of 
Grassmann variables and scalars in EqnJHUs not problematic, since we shall use this discussion 
solely to define a partition function for the entire system. And soon thereafter we shall 
integrate over the electron degrees of freedom, so that only a functional involving the scalar 
potential survives. 

The total Hamiltonian TC for the system, including the relativistic, quantum-mechanical 
part for the electrons is: 



TC = TLq + TCq 
H q =-Jl Q ft(-i^- V' + a)j> - fift $ (9) 

where 7 denote the Dirac gamma matrices, corresponding to the x, y, z directions. /1 is the 
electronic chemical potential, which shall be determined by demanding charge neutrality. 
The relativistic quantum correction has been encapsulated in the following dimensionless 
parameter: 

TZq = (10) 
Q Ak B T v ; 

where c is the speed of light. 

The parameter a in Eqn. M is related to the rest mass of the electron: 

mcA . . 

where m is the mass of the electron. It turns out that the chemical potential fi » k^T 
under conditions representative of white dwarfs. It follows that the relativistic approximation 
employed in Eqnj9] is valid. 



The partition function may be written in scaled variables as: 
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Z = J V4>V 2 ^ ± exp(- J d 3 x'(n + n Q )) (12) 

where it is understood that ksT = 1 in the units we are using. The quadratic nature 
of the energy functional in Eqn{12] allows us to perform the functional integration ove r 
the Grassmann variables associated with the electronic degrees of freedom (jRamondl Il98ll ) . 
allowing us to obtain: 



Z ~ J exp(- J dV((l/2)| V0| 2 + (l/2)0 2 - 0)) y/Det(F) 

Det(^) =exp(Trln(J 7 )) 

T = K Q (-Z7 • V + a) + AttT'4) - 0n (13) 

where (3 = 1/(/c_bT). Having integrated over the electronic degrees of freedom, we are left 
with an effective energy density in terms of the electrostatic potential alone. We could have 
also done things the other way, integrating over the electrostatic potential in the partition 
function, leaving a quartic in the fermionic variables, as is conventionally done. Our proce- 
dure can be said to have bosonized our plasma, since we now only have the scalar potential to 
investigate. We shall show below how our method leads to useful insights into the statistics 
of the plasma. 

We now need to evaluate the determinant of the operator obtained in the process of 
performing the quadratic functional integral over fermionic variables. This is conveniently 
performed in Fourier space, and to quadratic order in the electrostatic field variable: 



Tr ln(jF) ?a Tr J 



d 3 k 




(2tt) 3 

7 -k + a- Ph/Hq yy-k + a- /3fi/TZ Q 

where 4>(k) is the Fourier transform of and b is defined as: 

, 47rr' , . 

b=^~ (15) 

For temperatures in the range of 2 x 10 8 i^, and using a number density for our plasma 
~ 2 x 10 32 cm -3 (which corresponds to a mass density of 3 x 10 9 g — cm" 3 ), we find with 
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Z ~ 10, that A ~ 1 Fermz(10 _12 cm), T' ~ 10, and 7?.q ~ 10 3 . This is an indication that 
quantum corrections are dominant in this system. Hence AttT'/TZq ~ xl0~ 2 << 1 is an 
excellent choice for a perturbation parameter, so that the quadratic expansion considered in 
Eqnfj3]will be sufficient for our purpose. 

Terms devoid of field variables in Eqn. [TH will be ignored. For the particular case 
of the UDP considered here, the cubic term ignored in the above expansion is of O(10~ 2 ) 
compared to the bilinear terms we retained. Thus, higher order terms can be safely ignored. 
We must also consider for consistency higher order terms in the expansion of the Boltzmann 
distribution corresponding to the nuclear charges. We shall do so in the next section. 

We can now determine the chemical potential by demanding charge neutrality for the 
system. This is accomplished by considering the k — > limit of the terms in the energy 
functional that are linear in the electrostatic potential 0(fc), and setting the coefficient of 
the linear term to zero. 

= (8nV + mc 2 /k B T) » 1 (16) 

k B T 

Then the corresponding Lagrangian density in real space may be written to order b 2 as 
follows: 



^-effective = -dV0(r)| 2 



ffective 

(0(r)) 



V effecU ve(m) « 26 2 0(f)((V 2 - 4b 2 )-' - 8b 2 (V 2 - Ab 2 Y 2 )<P(r) (17) 

where b = 4^1"/ Aq. The equation of motion obtained by extremizing the above Lagrangian 
with respect to variations in the field \ is linear, albeit of order six: 



((A-46 2 ) 2 (l + A) + (A + 46 2 ))0 =0 

A EE V 2 

Upon factorizing, the trinomial in Eqn. [TS] may be cast as: 



(A-s 3 )(A-s 2 )(A- Sl )0(r) = O 



(19) 



- 8- 



where si,S2,S3 are the roots of the trinomial in Eqnfl8l These can be found easily using 
Mathematica. Upon careful examination of these roots, only one yields the correct limit for 
the screening length given by l/y/s v l/y/s 2 , l/y/s 3 in the limit that % — > oo. The screening 
length must go to infinity in this limit due to quantum fluctuations.- On e must remember 
that quantum fluctuations increase the screening length tjChitanvisI 120071 ). We will denote 
this root by s\. It may be written down simply as follows: 



si « 4.936 2 « 1 (20) 

For this limiting case, the dimensionless screening length is then \/4.93/6 >> 1. Hence 
we shall choose the physically interesting solution to satisfy: 

(V 2 - Sl )0(f) = (21) 

and this will be sufficient to guarantee that the sixth order differential equation is auto- 
matically satisfied. One may verify the screening length argument by noting that in one 
dimension exp(— 0:^4.93/6) solves the differential equation. We have used here a finite- 
temperature formalism involving not only a relativistic quantum electron gas, but also an 
admixture of positive nuclei to study screening in the limit that electronic quantum fluctu- 
ations dominate. In our case, the screening length turns out to be: 



^screening 



a 



A(T) 
2 1 



e 
he 



137 



(22) 



This screening length will be discussed in much detail in section 4. It may be compared 
to t he relativistic Thomas- Fermi s creening length at zero temperature for a pure electron 
gas Jjanovicilll977l : lltoh et alil2002h : 



A 



Thomas— Fermi 




5.2 x 10™(Zp 6 /Ay/3 
where p 6 is the mass density in powers of 10 6 , and A is the atomic weight. 



(23) 
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Upon inserting values for variables that are representative of white dwarfs, it turns out 
that Xscreening ~ 10~ 12 cm, while \ T homas- Fermi / ^screening ~ 10. We ascribe this difference 
to two factors, viz., finite temperature effects, and the fact that we are considering here a 
mixture of positively charged nuclei and electrons. One can see then that our approach will 
provide a higher fusion enhancement factor. 



3. Non-linear screening effects of nuclear charges 

In the previous section, we discussed in detail how non-linear terms arising from electron 
degeneracy are safely of much smaller magnitude than those retained. The issue appears to 
get turned around when one considers non-linear terms arising from higher order terms in the 
expansion of the Boltzmann distribution corresponding to nuclear charges. This is because 
the coefficients of the nonlinear terms are larger than those found for degenerate electrons 
in the previous section. The issue may be compactly discussed by noting that retention of 
the quadratic term in Eqn. [3] leads to the following cubic modification of EqnJHJ 

v -> v - -</> 3 (24) 
6 

The corresponding partition fu nction for our classical system resembles one for a Eu- 



clidean scalar quantum field theory (jRamondl 1 1 98 ll ) . Effects of the non-linear terms can be 



evaluated perturbatively, using Feynman diagrams. Feynman diagrams can be used to esti- 
mate the leading order contribution from this cubic term to the self-energy of the system (or, 
equivalently, the dielectric constant, or the screening length). It turns out that the actual 
value of the contribution is numerically quite small. This suggests that the linear screening 
approximation we retained in the previous section is an acceptable approximation. 

In order to perform this calculation, we formally ignored terms of 0(b) ~ 10~ 2 << (1/6) 
in our diagrammatics. We then used Mathematica to formulate symbolically the lowest order 
contributions from the cubic potential. The term that survives is a polarization-like diagram 
which comes from terms of (9((l/6) 2 ). The calculation is done as usual in momentum space. 
At this point the momentum variable is simply set to zero, so that we get an expression for 
the inverse square of the screening length: 



-* 1 f dp -» -» 

^{k)polarizatian-like{k) = — J Go(k)Go(p — k) 

Go(p) = (25) 
p z + 1 
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where is the usual self-energy. It can be easily related to the dielectric constant of the 
UDP. We shall restrict attention to k — 0, when the self energy it reduces to the inverse 
square of the screening length, and is sufficient to allow us to gauge its magnitude relative 
to the degenerate contribution in the previous section. 

Whence: 

^-polarization-like = S(& = 0) polarization-like = 

~ 3.3 x 10 -3 (white dwarfs) 
^poiarization-uke ~ 1.7 x 10~ n cm (white dwarfs) (26) 



We see that the screening length from the cubic term in the energy functional (due to 
nuclear charges alone) is larger by an order of magnitude than that obtained in the previous 
section by an order of magnitude (A screening ~ 10~ 12 cm for a white dwarf, from EqnJT9l). This 
difference in the screening due to nuclear and electronic charges means that electrons will be 
predominantly in closer proximity to a nucleus than another another nucleus. Since we are 
interested in screening at extremely short distances, we can therefore ignore the contribution 
to screening from Eqn. [26j 

The conclusions of this section may have to be modified if higher order terms turn out to 
be larger. For that purpose, we retain the quartic term in the energy functional arising from 
the Boltzmann distribution for the nuclear charges, so that Eqnj24l is modified as follows: 



(27) 



The lowest order contribution from the quartic term is the setting-sun diagram (JRamond 



198l|) . There are no cross-terms at this order between the cubic and quartic terms. 



1 f f d 3 h d 3 k 2 



>2_ 

6 J J (27r) 3 (27r) 3 ' 



^setUng-sun(p) = ~ / / , n _, a , n _, a Go(fcl)G (fc 2 )G (lp - h - k 2 \) (28) 



Using t^ 1 = J dXexp(-Xt), converting the momentum integrals to center of mass and 
relative co-ordinates, and using dimensional regularization, we can compute the self-energy 
for p = in the following form: 
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^setting— sun (0) 



1\ L I 



QJ 8tt 4 
°° exp(-fo) 



ln(£) - ln(l + y) , , , 

— ^ \J dy 29 

where £ — > + , a small-distance cut-off, and e — > + have been inserted to guarantee conver- 
gence. Using Mathematica, it can be shown that: 



L = ({)-T" 

/ = h^) _ 1 / 3- 4 cosec^) + 1; 2 _ £ , _ 



4e 16 

where 2-^1(0, b; c; z) is the hypergeometric function. 

Employing Laurent expansions within Mathematica in powers of e, and upon using 
counter-terms that account for divergences, one obtains the following correction to the square 
of the screening length, correct to lowest order in the quartic term: 



^setting— sun (0) 

^ (-7 - MO + 1 (-2^ + 12r' log (I)) - 1 ((2 V(3) - 7,))) 

V? = I (3D 

where a « A is a microscopic length cut-off required to render the integrals finite, V w 0.577 
is Euler's constant, and 71 ~ 0.572 is the value of the derivative of the Hypergeometric 
function which appears in Eqnf30l A logarithmic dependence of our answer on this cut-off 
implies a relative insensitivity to this parameter. It is clear that the theory used in this 
section is certainly not valid at the nuclear level, and so we will use £ = a/A ~ 10 2 . A 
cursory examination shows that this particular diagram yields a small contribution O(10 -4 ) 
contribution to the self-energy. Then following the argument above, for the cubic term, this 
screening contribution can be ignored as well. Thus there appears to be a trend for the 
classical, nonlinear terms to be small. 

Of course, further issues may arise in this perturbative argument as even higher order 
terms arising from the Boltzmann distribution are contemplated. We will leave these ques- 
tions for future investigation. It may be possible to extend to white dwarfs the methods 



utilized by Brown et al (lBrown et al 
highly charged plasma. 



20061 ) for obtaining the screening length for a dilute, 



4. Quantum effects of nuclei on screening 



Over the years, researchers have delved into the importance of applying a quantum- 
statistical treatment to the nuclei surrounding the ones undergoing fusion. This is a reason- 
able point to investigate, given the extremely high densities available in white dwarfs. Path 
integral Monte Carlo (PIMC) techniques have led to the discove ry that the bosonic nature 
of nucle hlltoh et al.lll990l ; IOgatalll997l ; iPollock and Militzerl 120041 ) makes a calculable, small 
contribution to the screening length, or equivalently, to the dielectric constant of the UDP 
under consideration. Here we provide an analytical underpinning to that observation. The 
second point that needs to be reinforced is that quantum effects of nuclei in the UDP are 
small in general, compared to the electronic contribution. The argument is basically that 
of Born and Oppenheimer, who showed that the nuclear mass is so large compared to the 
mass of an electron that an adiabatic approximation can be applied. That argument has to 
be extended to finite temperatures. 

We begin by considering a UDP consisting of spin-zero nuclei, in addition to a sea of 
degenerate electrons. Thus the Poisson-Boltzmann equation (EqnJH]) will be replaced by: 



v(4>, $) -> v B {4>, *, $±) = -47rZ0r'(^ t #) + 4vr0r'(> t ^) (32) 

where 4f', \fr are the second-quantized creation and annihilation operators corresponding to 
bosonic nuclei, and we have continued to assume that the photons are numerous and hot 
that they can be treated classically. 

Equation [9] must be modified to account for the free Hamiltonian of the nuclei: 



7~£q ^7~£q + T~(-AfQ 
H MQ = A Q (M) |W| 2 
/ h 2 A~ 2 \ 

A Q (M) = — — (33) 
QK ' \2Mk B T J v ; 

where Aq(M) is defined in terms of the mass M of the nucleus under consideration. It has 
been assumed that the interaction between nuclei are dominated by the Coulomb potential, 
and we have employed a non-relativistic treatment for nuclei. 
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The corresponding functional integrals involving the nuclear field variables can be done, 
just as the fermionic degrees of freedom were accounted for. The net result for the partition 
function is: 



where 



Z~ Jv4> exp(- J (fx (1/2)|V0| 2 ) y/Det(F) Det(S)" 1 



(34) 



Tr ln(£ _1 ) 



d 3 k 
(2tt)3 



ln(4vrZr'0(A;) + A Q (M) k 2 



(35) 



Using parameters representative of a white dwarf, and given that M ~ 2 x 10 4 m, Z ~ 10, 
it turns out that A Q (M) ~ 10" 4 A Q (m), so that AQ(M)/4vrZr' ~ 10~ 3 . This clearly 
shows that quantum effects, including Bose statistics of nuclei in the UDP are negligible 
in circumstances representative of a white dwarf. Thus, we may continue to use a classical 
treatment for nuclei, as was done in the previous section. One can regain from Eqn l35l the 
Boltzmann approximation using the method outlined by Chitanvis rtChitanvisI 120071 ). This 
will yield systematic, extremely small quantum corrections to the Boltzmann approximation 
in powers of h 2 . Considering the nuclei as fermions does not change this conclusion. 

Quantum effects of nuclei are much smaller than the quantum effects due to electrons, 
primarily due to the large mass difference. As such they can be ignored. Our c onclusions 



conce rning th e effects of i ndistinguishability are in general agreement with Itq h fjltoh et al. 



19901 ). Ogata (10gatalll997r ) and Pollock and Militzer flPollock and Militzerll2004h . These au- 
thors did not consider the effects of electron degeneracy on the same footing as the nuclear 
quantum effects. 



5. Results 



We shall compare our results with those found in published literature. For succinctness, 
we shall quote the enhancement f actors of pycnonuclear reactio n rates given in the review 
ar ticle of Ichimaru and Kitamura tllchimaru and Kitamural 119991 ) and results by Gasques et 
a KlGasques et al.l 120051 ) . 



Our approach provides an integrated, first-principles, systematic theory of screening 
effects due to electrons and associated nuclei. As such, it is possible for us to estimate the 
accuracy with which we can calculate our enhancement factors. In order to do that for 
white dwarfs, we must first generalize our formalism to a binary ionic mixture (BIM), e.g. 
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a mixture of 12 C, 16 nuclei and associated electrons, as representative of the species in a 
white dwarf. This is easily accomplished via the following substitutions in all our formulae: 



Z 2 _^ 2 2 = niZ ^ + U2Z * 

n 

n — > n = n\ + ri2 (36) 

where ni, n 2 are the average number densities of each species in the UDP. These substitutions 
arise naturally through a re-derivation of our theory for a BIM. Generalizations to more than 
two components is straightforward. 

Ichimaru uses improvements over a standard procedure to obtain net enhancement fac- 
tors. First, the screening length/dielectric constant of an electron gas in a jellium of positive 
ions is obtained. Then various sophisticated methods are used separately to obtain screening 
effects due to nuclei surrounding the moieties undergoing fusion. Physically reasonable mix- 
ture rules are utilized to obtain the overall enhancement of nuclear rates in a UDP, caused by 
the screening of the nuclear Coulomb repulsion by intervening charges. It is possible to gauge 
the accuracy of such calculations by using our systematic approach, where the magnitude of 
terms neglected can be estimated. In general we find our estimates for fusion rate enhance 



ment in white dwarfs do agree with those of Ichimaru and Kitamura fllchimaru and Kitamura 



19991 ). One case where there is severe disagreement can be attributed to the breakdown of 
approximations made in this paper (see Table [T|). The values of the classical plasma pa- 
rameter Tdassicai = Z 2 e 2 /r Q kT « 170 (r is the mean-free distance between particles) for 
the cases listed in Table [H As such there is no concern regarding the UDP being close to a 
crystallized state. 



Gasques et a KIGasques et al.ll2005f) utilize a re-par ameterized version of the enhancement 



factor obtained by Slattery et a KjSlattery et al.lll982l ). The enhancement factor of Gasques 



et al is defined in terms of the plasma parameter of a one-component classical plasma. In 
our case, we have a binary mixture. It is nevertheless instructive to compare to our results 
the classical enhancement using a positively charged gas, with each nuclei carrying a charge 
Ze. In general, there is a reasonable agreement. 

We note in passing that in order to obtain agreement between different approaches 
for calculating the fusion enhancement factor, it was essential to assume a relativistic de- 
scription for the degenerate electrons in a white dwarf, while simultaneously considering a 
background of positive charges representing nuclei, at finite temperature. Assuming a non- 
relativistic description for electrons caused us to obtain enhancement factors that are much 
more conservative. 
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The numerical values obtained for our screen ing length and comparisons to Ich imaru's 



enha ncement factor (jlchimaru and Kitamuralll999l ) and that of a classical gas of nuclei (IGasques et al 



20051 ) have been encapsulated for white dwarfs near ignition in Table [TJ 
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Quantity 


Current theory 


(Ichimaru and Kitamura 1999) 


(Gasques (2005)) 


Model 


WD1 






Density 
T (Kelvin) 
Composition 


<- 
<- 


3.0 x 10 9 g - cm" 3 
1.8 x 10 8 K 
12 C, e~ 


-> 
-> 


^enhance 


9.02 


12.09 


11.84 


Model 


WD 2 






Density 
T (Kelvin) 
Composition 


<- 
<- 


9.0 x 10 9 g — cm' 3 

1.1 x 10 8 ir 
12 C, e" 


-> 
-> 


1 enhance 


32.7 


23.10 


27.96 


Model 


WD 3 






Density 
T (Kelvin) 
Composition 


<- 


9.0 x 10 9 g — cm -3 
3.4 x 10 7 K 

12 C, e~ 


-> 


^enhance 


90.45 


20.76 


190.3 


Model 


WD4 






Density 
T (Kelvin) 
Composition 




9.0 x 10 9 g — cm -3 

1.1 x 10 8 K 

12 C (75%), 16 (25%)e" 




^enhance 


36.60 


23.12 


31.62 



Table 1: Comparison of corrected rate fact ors for white dwarfs near ignition, between 



our c alculation and two previous calculations ( jlchimaru and Kitamura! 11999c iGasques et al. 
20051 ). The enhancement of the fusion rate is calculated as exp(r enhance), ^enhance = 
e 2 / (A screening k B T). Our quantum-influenced screening length A screening is defined in dimen 



sionless terms as l/^fs 1 via EqnJT9l The differe nt scenarios for white dwarfs , viz., models 
WD1 — WD4 are described in Ichimaru's paper (jlchimaru and Kitamuralll999l ). The reason 
for the major discrepancy for model 3 is that the temperature for this case is sufficiently low 
so that the approximations made in this paper do not hold. 



